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Abstract:  WIPL-D  is  being  used  to  compute  monostatic  and  bistatic  radar  cross  sections  of  a  trihedral  corner  reflector 
over  the  frequency  range  of  1-12  GHz.  Initial  results  are  discussed  for  vertically  and  horizontally  polarized  fields  at  two 
frequencies  (1.5  and  3.8  GHz).  These  computations  are  being  used  to  provide  a  benchmark  against  which  the  performance 
of  WIPL-DP,  a  parallelized  version  of  the  WIPL-D,  may  be  compared.  Recently,  sea-scatter  data  was  collected  from  1.9- 
11.5  GHz  with  trihedral  reflectors.  Since  the  WIPL-D  calculations  discussed  here,  as  well  as  planned  calculations,  will  be 
used  to  provide  predictions  of  sea  scatter  that  will  be  compared  to  the  actual  data,  the  sea-scatter  geometry  is  described. 
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1.  Introduction 

Monostatic  and  bistatic  radar  cross  sections  of  a  trihedral  corner  reflector  are  calculated  with  the  WIPL-D 
electromagnetics  code  at  two  frequencies  (1.5  GHz  and  3.8  GHz)  for  vertically  and  horizontally  polarized  fields 
to  provide  a  benchmark  against  which  WIPL-DP,  a  parallelized  version  of  the  WIPL-D  [1],  may  be  compared. 
WIPL-DP  is  currently  being  developed  under  the  sponsorship  of  the  DoD’s  High  Performance  Computer 
Modernization  Office  (HPCMO  -  http://www.hpcmo.hpc.mil/)  in  one  of  its  Common  HPC  Software  Support 
Initiatives  (CHSSIs),  entitled  Parallel  Scene  Generation/Electromagnetic  Modeling  of  Complex  Targets  in 
Complex  Clutter  and  Propagation  Environments.  This  CHSSI  is  a  joint  university-industry-govemment  effort 
led  by  the  Air  Force  Research  Laboratory  and  Black  River  Systems  Corp.  The  motivation  behind  the  trihedral 
analysis  is  to  provide  one  of  two  sufficiently  difficult,  but  pertinent  and  numerically  accessible,  Navy 
applications  to  test  the  capability  of  WIPL-DP.  The  relevant  Navy  problem  associated  with  the  trihedral  is  to 
understand  the  target-like  artifacts  in  radar  returns  called  sea  spikes  that  are  induced  by  ocean  scatter  [2]. 

Historically,  the  trihedral  comer  reflector  has  been  used  for  instrumentation  purposes,  radar  calibration,  and 
devices  with  enhanced  radar  cross  section  (RCS)  [3-5].  Until  recently,  most  scattering  studies  of  the  trihedral 
have  been  understood  through  optics  and  empirical  evidence,  as  a  result  of  the  excessively  large  computational 
burden  imposed  by  this  problem.  With  advances  in  simulation  tools  and  computational  capabilities,  the  trihedral 
comer  reflector  is  revisited. 

Trihedral  comer  reflectors  have  recently  been  used  in  a  series  of  ultrawideband  (UWB)  measurements  at  the 
Atlantic  Underwater  Test  and  Experiment  Center  (AUTEC)  in  the  Bahamas  to  investigate  the  low-elevation 
(grazing  angles  less  than  4°)  scatter  of  RF  signals  from  the  open  ocean  in  an  attempt  to  understand  and  mitigate 
the  phenomenon  known  as  sea  spikes  [6].  Sea  spikes  introduce  undesirable,  false,  target-like  artifacts  in  the 
received  signals  of  ship-based  radars  -  a  very  serious  impediment  to  detecting  real  targets.  A  state-of-the-art, 
UWB,  ultra-high  resolution,  dual-polarized  radar,  designed  and  built  by  Pete  Hansen  of  the  Radar  Division  of 
the  Naval  Research  Laboratory,  has  been  taking  these  short-pulse  UWB  measurements  at  S-band  and  X-band. 

The  AUTEC  geometry  (Fig.l)  has  a  1-m,  anodized-aluminum,  right-trihedral  reflector  mounted  on  a  tower 
that  is  8.1  m  above  the  sea  surface  and  242  m  from  the  shore-mounted  UWB  radar,  called  the  microwave 
microscope  (MWM).  The  transmitted  signal  travels  on  four  paths  before  it  returns  to  the  receive  antenna:  (1)  the 
direct  path  between  the  MWM  and  the  trihedral;  (2)  the  one-bounce  path  from  transmitter  to  ocean  to  trihedral 
to  receiver;  (3)  the  one-bounce  path  from  transmitter  to  trihedral  to  ocean  to  receiver;  and  (4)  the  double-bounce 
path  from  transmitter  to  ocean  to  trihedral  to  ocean  to  receiver.  For  this  work,  the  sea  surface  is  assumed  to  be  a 
perfectly  conducting  flat  plate,  so  that  one  need  only  consider  the  specular  reflections  corresponding  to  these 
four  paths.  Paths  (1)  and  (4)  are  monostatic  from  the  target’s  perspective,  while  paths  (2)  and  (3)  are  bistatic 
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relative  to  the  target.  Conventional  optical  analysis  has  only  been  applied  to  the  monostatic  radar  cross  section 
(RCS),  and  most  bistatic  RCS  has  been  obtained  by  experiment.  By  leveraging  the  computational  capabilities  of 
modem  workstations,  monostatic  and  bistatic  radar  cross  sections  (RCSs)  of  the  trihedral  reflector  are 
simulated. 

The  MWM  generates  an  UWB  short-pulse  signal  by  impulse  exciting  a  TWT.  Three  TWTs  are  used  to 
cover  the  frequency  bands:  1-2  GHz;  1.9-5. 5  GHz;  6.5-11.5  GHz.  The  nominal  center  frequencies  and  range 
resolutions  corresponding  to  these  bands  respectively  are  (1.5  GHz,  15  cm),  (3.8  GHz,  3  cm),  and  (9  GHz,  2 
cm).  Although  calculations  are  being  made  for  each  frequency  band,  only  calculations  at  1.5  GHz  and  one 
calculation  at  3.8  GHz  were  available  when  preparing  this  paper. 


Trihedral 


2.  Trihedral  and  WIPL-D  Geometries 

The  trihedral  is  composed  of  three  right  isosceles  triangles  (Anin3n8,  An3n2n3,  and  Aninxn2  in  Fig.  2),  where 
Fig.  2  shows  the  screen  capture  of  the  WIPL-D  interface,  the  largest  edges  (n2n3,  n3n8,  n8n2)  have  length  a 
(which  is  1  m  in  the  ensuing  calculations),  {n,}  are  the  trihedral  grid  points  from  WIPL-D,  and  P(n,)  is  the 
coordinate  of  n;.  WIPL-D  and  WIPL-DP  model  plate  structures  with  quadrilateral  surfaces  and  impose  the 
restriction  that  the  quadrilaterals  be  convex.  Therefore,  each  triangular  side  is  modeled  using  three  smaller 
plates.  To  ease  the  geometric-modeling  complexity,  the  three  shared  edges  (nin2,  nin3,  and  nin8)  between 
adjacent  triangular  sides  are  aligned  respectively  with  the  Cartesian  axes  x,  y,  and  z,  where  the  origin  is  placed 
at  vertex  n,  of  the  trihedral.  Spherical  angles  in  the  WIPL-D  codes  are  measured  from  the  x-axis  towards  the  y- 
axis  (azimuth  or  longitude  (f)  and  from  the  xv- plane  towards  the  z-axis  (elevation  or  latitude  6),  which  is  90°  - 
9,  where  the  co-latitude  9  is  the  usual  spherical  angle.  In  Fig.  2,  one  is  looking  into  the  trihedral  along  the  axis 
of  symmetry  through  vertex  ni.  In  3-dimensional  space,  ni  lies  behind  the  plane  of  triangle  Anin2n8 ,  and  the  axis 
of  symmetry  is  perpendicular  to  this  plane. 

All  three  triangles  that  comprise  the  trihedral  are  congruent  and,  consequently,  have  the  same  dimensions. 
In  particular,  from  a  planar  view  of  triangle  An,n2n3  (Fig.  3),  one  can  see  that  edge  n2n3  determines  the 
dimensions  of  Anin2n3,  edges  nin2  and  n,n3  have  length  0.707a,  point  n4  is  the  centroid  of  An,n2n3,  and  points  n5, 
n6,  and  n7  are  the  midpoints  of  their  respective  edges. 
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p(ni)  =  (0,0,0) 

P(n2)  =  (fl/V2,0,0) 

P(n3)  =  (0,a/V2,0) 

P  (  n4 )  =  (  «/( 3  V2 ),  a/(3  V2 ),  0  ) 
P(ns)  =  (0,a/(2>/2),0) 
P(n6)  =  (a/(2,/2),0,0) 
P(n7)  =  (a/(2V2),a/(2V2),0) 
P(n8)  =  (0,0,a/>/2) 

P(n9)  =  (0,0,a/(2V2)) 
P(n10)  =  (a/(3V2),0,a/(3V2)) 
P(nn)  =  (a/(2>/2),0,a/(2>/2)) 
P(n12)  =  (0,a/(3V2),a/(3V2)) 
P(n13)  =  (0,a/(2V2),a/(2V2)) 


Fig.  2.  Frontal  view  of  a  right  trihedral  corner  reflector  along  the  axis  of  symmetry  through  vertex 
nj,  with  the  complete  WIPL-D  node  geometry  and  scatterer  coordinate  system.  The  positive  .r,  y,  and 
z  axes  coincide  with  line  segments  nin2,  nin3,  and  n,n8,  respectively,  and  Pin,)  is  the  coordinate  of  n,. 


Fig.  3.  Geometry  of  the  single  triangular  side  An |ii2n3  comprising  the  trihedral,  where  P(n,)  is  the 
coordinate  of  n„  the  positive  x  and y  axes  respectively  coincide  with  line  segments  niii2  and  nin3,  and 
the  positive  z  axis  points  out  of  the  paper. 
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To  relate  the  AUTEC  and  WIPL-D  geometries,  consider  the  side  view  of  the  trihedral  (Fig.  4)  relative  to 
the  AUTEC  geometry  (Fig.  1).  In  this  configuration,  Of is  the  angle  from  the  vv- plane  to  a  horizontal  plane  ( 0 — 
Of)  parallel  to  the  ideal  sea  surface.  Note  that  this  plane  includes  the  trihedral’s  axis  of  symmetry  {(j)-  45°,  0—  Of 

-  35.26°).  The  angles  ft  and  a  are  specific  to  the  AUTEC  configuration,  where  ft  is  the  angle  between  the 
symmetry  axis  and  the  radar-trihedral  line  of  sight  {(j)  -  45°,  0=  Of  +  ft)  and  a  is  the  angle  between  the  radar- 
trihedral  line  of  sight  and  the  specular  point  on  the  sea  surface  (. <j>-  45°,  0-  0f  +  ft  +  a).  From  trigonometry,  a 

-  2.03°,  ft-  0.90°,  and  the  grazing  angle  is  a  +  ft-  2.93°. 


3.  WIPL-D  Experiments 

Two  basic  classes  of  numerical  experiments  have  been  performed:  monostatic  RCSs  at  1.5  GEIz  and  3.8 
GEIz;  and  bistatic  RCSs  at  1.5  GEIz  at  specific  incidence  and  arbitrary  observation  angles.  These  two  classes  are 
studied  in  anticipation  of  numerical  experiments  that  are  specific  to  collected  data  with  the  AUTEC  geometry 
and  future  RCS  measurements  of  trihedrals  in  a  compact  range.  Since  the  AUTEC  data  has  been  collected  at 
horizontal  and  vertical  polarizations  from  1-11.5  GEIz,  five  frequencies  of  interest  for  numerical  simulation  are 
listed  in  Table  1,  with  the  associated  WIPL-D  modeling  requirements  for  a  1-m  comer  reflector. 

Despite  its  simple  geometric  structure,  the  trihedral  comer  reflector  is  a  complex  target  for  computational- 
electromagnetics  solvers.  As  the  frequency  and  target  size  a  increase,  the  number  of  unknowns  required  by 
WIPL-D  grows  exponentially.  This  scaling  in  computational  requirements  quickly  poses  problems,  even  for 
modem  workstations,  and  indicates  that  the  trihedral  comer  reflector  is  an  excellent  benchmark  for 
computational  solvers.  The  storage  requirements  are  independent  of  the  computer  used;  however,  solution  times 
will  vary  with  computer.  The  results  shown  here  are  computed  on  a  Pentium-4,  2.5-GHz  workstation  with  1  GB 
of  RAM  and  512  kB  of  L2  cache.  Unfortunately,  the  storage  required  to  simulate  the  7-1 1  GHz  range  is  beyond 
the  capability  of  any  32-bit  workstation.  Thus,  to  collect  data  at  these  higher  frequencies,  a  more  sophisticated 
computer  must  be  used  -  for  example,  a  parallel  computer  or  64-bit  workstation.  Consequently,  only  data  at  1.5 
and  3.8  GHz  can  be  shown  at  this  time. 
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Frequency 

all 

Number  of  Unknowns  (u) 

Storage  Required  (8u2) 

1.5  GHz 

5.00 

960 

7.2  MB 

3.8  GHz 

12.67 

5131 

200  MB 

7.0  GHz 

23.33 

13955 

1.5  GB 

9.0  GHz 

30.00 

25702 

5.1  GB 

11.0  GHz 

36.67 

35833 

9.8  GB 

Table  1.  WIPL-D  modeling  requirements  for  various  frequencies. 


3.1.  Monostatic  RCS 

With  the  trihedral  aligned  along  the  x,  y,  and  z  axes,  the  trihedral  acts  as  a  comer  in  the  octant  containing 
positive  x,  y,  and  z  components,  which  corresponds  to  0°  <  (f>  <  90°  and  0°  <  6  <  90°.  While  most  of  the 
interesting  scattering  occurs  in  this  octant  (the  optical  field  of  view  of  the  trihedral),  some  scattering  does  occur 
in  the  other  seven  octants.  In  Fig.  5,  three  scattering  zones  are  apparent  in  the  resulting  3-dimensional  RCS 
patterns  of  the  vertically  polarized  field  (horizontal  and  vertical  refer  to  the  ^-component  and  the  ^-component 
of  the  field,  respectively).  This  pattern  is  obtained  by  mapping  the  nonnegative  RCS  value  as  the  distance  from 
origin  ni  in  the  (<p/J)  direction.  The  strong  region  in  the  middle  of  the  pattern  occurs  when  incident  waves  are 
reflected  off  all  three  triangular  faces  -  hence  this  area  is  termed  the  triple-bounce  zone.  Similar  scattering 
regions  occur  when  waves  bounce  off  two  sides  ( double  bounce )  and  one  side  {single  bounce).  The  peak  RCSs 
are  seen  in  the  triple-bounce  zone,  which  is  the  dark  hexagonal  region  centered  about  the  trihedral  axis  of 
symmetry.  Beyond  the  hexagonal  region,  the  RCS  drops  off  quickly,  until  it  significantly  increases  for 
directions  that  are  nearly  parallel  to  the  trihedral  boundaries.  The  triple-bounce,  hexagonal  region  is  called  the 
effective  aperture  [5]. 


+y 


+Z  +2 

Fig.  5.  Three  views  of  the  3-dimensional  vertically  polarized  monostatic  RCS  for  a  1-m  trihedral  at 
1.5  GHz  within  the  optical  field  of  view  of  the  positive  xyz  octant,  where  the  data  are  normalized 
relative  to  l2  and  the  dark  and  light  shading  denote  large  and  small  RCS,  respectively. 


double 

bounce 


single 

bounce 


bounce 


The  monostatic  RCS  at  1.5  GHz  is  plotted  in  Fig.  6  for  two  perpendicular  planes  that  intersect  the  surface  in 
Fig.  5:  (a)  45°  and  0°  <  0<  90°;  and  (b)  9-  35.26°  and  0°  <  <p<  90°.  In  Fig.  6b,  the  vertical  and  horizontal 
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polarizations  are  symmetric  about  (j)  =  45°,  as  expected  from  the  geometry,  and  decrease  slowly  as  <f>  moves 
away  from  45°  for  25°  <  (f>  <  75°,  where  they  differ  by  less  than  2  dBsm.  Near  <f>  -  12°  and  78°,  the  horizontal 
RCS  drops  significantly  (more  than  5  dBsm).  As  (j)  approaches  0°  and  90°,  both  polarizations  increase  to  a 
significant  fraction  of  the  peak  values  attained  along  the  axis  of  symmetry  {(j)  =  45 °,9  =  35.26°),  because 
significant  backscatter  is  produced  by  triangles  Anin8n2  and  Anin3n8,  respectively. 

For  case  (a),  where  the  plane  symmetrically  bisects  the  trihedral,  the  RCSs  for  both  polarizations  behave 
asymmetrically  as  9  varies  from  0°  and  90°  (Fig.  6a),  that  is,  as  9  varies  from  direction  nin7  to  direction  nin8. 
Again  the  polarizations  tend  to  track  each  other  to  within  3-5  dBsm.  In  addition,  both  polarizations  have 
significantly  lower  RCSs  near  9-  74°,  where  very  little  energy  (approximately  1/100)  is  reflected.  When  the 
plots  in  Fig.  6  are  visually  compared  to  results  in  the  literature  [4, 5, 7-9],  the  curves  have  the  same  general 
shape,  with  the  curves  in  Fig.  6  having  more  ripple  than  those  predicted  by  optical  theory  [7-9],  because  WIPL- 
D  includes  the  effects  of  diffraction.  For  example,  in  Figs.  9  and  10  of  Ref.  [8],  except  for  differences  in 
magnitude  of  the  local  minima  and  the  location  of  one  null,  the  RCSs  for  horizontal  polarization  appear  to  be 
smoother  versions  of  Figs.  6a  and  6b,  respectively,  for  comparable  a! A  (4.99). 


Fig.  6.  Monostatic  RCS  of  trihedral  at  1.5  GHz  for  vertically  polarized  (dotted  curves)  and 
horizontally  polarized  (solid  curves)  incident  fields  calculated  with  WIPL-D:  (a)  (/)  =  45°  and  0°  <  0 
<  90°;  (b)  9=  35.26°  and  0°  <  90°. 


In  Fig.7,  the  monostatic  RCSs  are  displayed  for  1.5  GHz  and  3.8  GHz  as  2-dimensional  projections  of  the 
3-dimensional  surface  mappings  (like  Fig.  5)  onto  the  plane  that  is  perpendicular  to  the  trihedral  axis  of 
symmetry.  The  positive  octant  (0°  <  </><  90°  and  0°  <  0<  90°)  is  mapped  onto  the  equilateral  curved  triangles  of 
Fig.  7,  where  white  and  black  of  the  gray-scale  bar  indicate  high  (20  dBsm)  and  low  (0  dBsm)  RCSs, 
respectively.  On  close  scrutiny  of  the  1.5 -GHz  plots,  one  can  discern  the  hexagonal  region  associated  with  the 
effective  aperture  (triple-bounce  region)  that  returns  the  greatest  energy  in  the  optical  approach  and  the  low 
RCS  returns  from  directions  near,  but  not  equal  to,  the  Cartesian  axes  ( x,y,z ).  Moreover,  the  horizontal  and 
vertical  RCSs  at  1.5  GHz  have  noticeable  differences,  for  example,  the  location  and  depth  of  the  nulls.  In 
contrast,  at  3.8  GHz,  the  effective  aperture  is  considerably  bigger,  with  slightly  larger  RCSs  for  the  horizontally 
polarized  fields,  and  the  difference  between  polarizations  is  less  obvious. 
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Fig.  7.  Two-dimensional  monostatic  RCS  plots  of  trihedral  calculated  with  WIPL-D:  (a)  horizontal 
polarization  at  1.5  GHz;  (b)  vertical  polarization  at  1.5  GHz;  (c)  horizontal  polarization  at  3.8  GHz; 
and  (d)  vertical  polarization  at  3.8  GHz. 


3.2  Bistatic  RCS 

For  the  bistatic  case,  RCS  is  plotted  versus  arbitrary  receive  (observation)  angles  (^0bs>#obs)  in  the  positive 
xyz  octant  for  three  incident  angles  along  the  geometric  half-plane  of  symmetry  {<(>  —  45°);  (1)  9-  20.26°,  which 
is  15°  above  axis  of  symmetry  (Figs.  8a, d);  (2)  9-  35.26°,  which  is  the  axis  of  symmetry  (Figs.  8b, e);  (c)  9  = 
50.26°,  which  is  15°  below  axis  of  symmetry  (Figs.  8c, f).  On  studying  Fig.  8,  one  may  infer  that  the  bistatic 
RCS  is  significantly  smaller  than  the  monostatic  RCS,  that  the  patterns  are  symmetric  about  the  geometric  plane 
of  symmetry,  and  that  the  polarizations  yield  different  patterns.  As  0obs  increases,  the  area  of  greatest  RCS 
travels  down  the  bisector  of  the  curved  triangular  region  and  transforms  from  being  vertically  elongated  to 
being  horizontally  elongated,  regardless  of  the  polarization.  For  0„\K  near  the  monostatic  direction  (35.26°),  this 
area  corresponds  to  the  trihedral’s  effective  aperture.  Also,  note  the  well  defined  bright  spots  outside  the  main 
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area  of  large  RCS,  which  indicate  that  one  can  expect  enhanced  RCS  from  regions  other  than  the  main  area. 
How  the  preceding  results  relate  to  the  bistatic  angle  is  left  for  a  future  discussion.  To  get  a  more  concrete  idea 
of  the  differences  between  polarizations,  three  cross-sectional  curves  are  plotted  for  three  incident  directions 
and  two  receive  angles  (0obs  =  35.26°,  ^obs  =  45°)  in  Fig.  9. 
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Fig.  8.  Two-dimensional  projections  of  bistatic  RCSs  of  1-m  trihedral  at  1.5  GHz  for  three  incident 
directions  (jfecAc)  and  both  polarizations  (horizontal  -  HPol,  vertical  -  VPol):  (a)  HPol  @ 
(45°, 20.26°);  (b)  HPol  @  (45°, 35.26°);  (c)  HPol  @  (45°, 50.26°);  (d)  VPol  @  (45°, 20.26°);  (e)  VPol 
@  (45°, 35.26°);  (f)  VPol  @  (45°, 50.26°). 


4.  Summary 

This  work  represents  initial  WIPL-D(P)  computational  modeling  of  monostatic  and  bistatic  scattering  from 
the  trihedral  comer  reflector  as  part  of  a  larger  application  problem  to  characterize  radar  sea  scatter.  It  is 
important  to  perform  additional  WIPL-D(P)  calculations  for  several  frequencies  from  1-12  GHz  so  that:  (1)  it 
may  be  determined  through  convergence  analyses  whether  WIPL-D(P)  produces  the  best  answer  it  can,  given 
the  limitations  of  the  computer  used  to  make  the  calculations;  and  (2)  WIPL-D(P)  predictions  can  be  compared 
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to  existing  UWB  radar  data.  In  addition,  the  trihedral  reflectors  used  for  the  collected  radar  data  need  to  be 
measured  in  a  compact  range,  and  these  range  measurements  should  be  compared  to  WIPL-D(P)  predictions. 


Fig.  9.  Bistatic  RCSs  of  1-m  trihedral  at  1.5  GHz  for  incident  directions  [(45°, 20.26°),  (45°, 35.26°), 
(45°, 50.26°)]:  (a)  HPol  @  6>obs  =  35.26°;  (b)  VPol  @  9oh 8  =  35.26°;  (c)  HPol  @  ^obs  =  45°;  (d)  VPol 
@  < boh 8  =  45°. 
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